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ABSTRACT 



A Ready Supply Store (RSS) containing repair parts which 
are anticipated to be used during the production process has 
been established to support the Naval Air Rework Facility 
(NARF) . While this supporting inventory has previously been 
constructed using historical demand data, a model which com- 
putes stock levels based on the next quarter's production 
schedule has been proposed. This thesis extends the use of 
projected production information in calculating RSS inventory 
levels from one to two periods, and compares the expected 
total costs from both systems under the assumption of a bi- 
nomial demand distribution which is appropriate to a NARF. 

As a result of this comparison, the conclusion is made that 
the two-period model offers only very minor expected cost 
advantages over a single-period formulation, while also being 
much more difficult to utilize due to the complex calcula- 
tions involved in the computation process. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Naval Air Rework Facilities (NARFs) perform overhaul and 
maintenance actions on various components utilized by fleet 
aviation units. These rework activities are accomplished in 
accordance with a quarterly production schedule which is 
established by a joint workload conference between the NARF 
and the Naval Air Logistics Command. Information on antici- 
pated NARF workload requirements is also available for sever- 
al additional future periods. Accordingly, the. NARFs are in 
the process of implementing a Material Requirements Planning 
(MRP) system which will utilize the available forecast to 
project requirements for individual spare parts used in the 
rework process. By establishing such a system, the NARF will 
be able to accomplish its assigned mission more efficiently 
by reducing tne number of work stoppages caused by stockouts 
[Ref. 1: 10]. 

MRP systems are intended to reduce or even eliminate re- 
pair part inventory requirements through the technique of 
phasing item arrivals to coincide with their need within the 
production process. However, variations in procurement lead- 
time and in the actual production schedules necessitate the 
establishment of some form of backup inventory support. For 
this purpose, the Naval Supply Center (NSC) supporting the 
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NARF has established a Ready Supply Store (RSS) of parts 
which are anticipated to be used during the production pro- 
cess. Once located in the RSS, these parts are protected 
from issue to other NSC customers. Consequently, factors 
such as the need to maintain maximum visibility of critical 
aviation parts within the entire Navy Supply System, the 
limitations on funding available to the RSS, and the require- 
ment to support customers other than the NARF mandate that 
the range and depth of items within the RSS inventory be 
carefully monitored. 

Prior to development of the MRP system, RSS stock quanti- 
ties were established by the traditional method of using his- 
torical demand data to compute appropriate levels for 
individual items. This technique is effective in situations 
where the workload is reasonably constant from one period to 
the next, but can result in large surplus or shortage quanti- 
ties when the production schedule changes to any great ex- 
tent. Since component overhaul activity at the NARF is 
susceptible to quarterly fluctuations, a system which deter- 
mines phased repair part requirements based on anticipated 
production schedules offers significant advantages over one 
founded solely on demand history. Accordingly, an appro- 
priate single-period inventory model utilizing the known pro- 
duction information for the next quarter was proposed by 
McMasters [Ref. 2: 4-15]. This model distinguishes between 

items which are replaced 100% of the time during overhaul and 
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those which are not, focusing on the latter category since 
the former should necessarily be stocked in quantities al- 
lowing for total replacement. The ensuing inventory system 
is then based on obtaining a balance between expected short- 
age and surplus costs in determining optimal stock quantities. 

B. PURPOSE 

While the McMasters single-period model represents a sig- 
nificant step forward by utilizing projected future require- 
ments in place of historical demand, it does not make maximum 
use of available production forecasts since only the workload 
for the next quarter is considered. In particular, surplus 
costs projected at the end of a given period on a specific 
item may be greatly reduced by considering the anticipated 
demand for that item in future periods. Alternately, a sig- 
nificant demand forecast in ensuing periods will probably 
result in a decreased tendency on the part of the system to 
accept a potential shortage of that item in the initial 
period. A model which incorporates anticipated rework acti- 
vity beyond the upcoming production period should provide a 
more accurate reflection of total expected costs over a par- 
ticular length of time and, thus, offers the potential of 
creating a more cost effective inventory mix. 

Accordingly, the objectives of this thesis are to: 

1. Develop a model which utilizes available production 
schedules for the next two quarters in determining 
appropriate inventory quantities for the supporting 
RSS. 
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2 . Conduct a cost comparison between the optimal results 
obtained from single and dual-period models. 

Completion of this two-period inventory system will repre- 
sent the first step toward development of a multi-period 
model which considers all available production information 
and covers a full procurement leadtime for each RSS stock 
item. 

C. PREVIEW 

Chapter II describes the general model utilized for the 
duration of this paper and its behavior when used in a 
single-period analysis. In Chapter III, this basic model is 
used as a two-period inventory system whose performance is 
examined under the assumption of a uniform demand distribu- 
tion. Chapter IV describes the two-period model under the 
assumption of a binomial demand distribution. This distribu- 
tion provides the best representation of demand appropriate 
to a NARF [Ref. 2: 4]. A sensitivity analysis on the re- 

sults obtained from using the two-period binomial model is 
conducted in Chapter V. Finally, Chapter VI provides a sum- 
mary of results as well as recommendations for future actions. 
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II. THE SINGLE-PfiRIOD MODEL 



A. PRELIMINARY INFORMATION 

An inventory model compatible with the general ideas dis- 
cussed in the first chapter must possess several essential 
features. These are: 

1. The model should select a stock quantity for each spe- 
cific item which minimizes inventory-related costs. 

2. As with the McMasters single-period system, the model 
should evaluate the trade-offs between expected surplus 
and shortage costs in computing appropriate inventory 
levels . 

3. The model should be able to incorporate two or more 
production periods in the planning horizon. 

4. The model should possess the capability of dealing with 
production schedules which vary from one period to the 
next. In other words, it must be able to treat demand 
probability distributions which have different parame- 
ter values in successive periods. 

The general system described in the remainder of this thesis 

possesses the essential features detailed above and is an 

adaptation of a model developed by Karlin [Ref. 3: 231-243]. 



B. GENERAL ASSUMPTIONS 

Before proceeding to examine the single-period model, it 
is necessary to specify the general assumptions which have 
been made in its formulation. These assumptions are listed 
as follows: 

1. No backorders are permitted. If the required spare is 
not available from RSS stock when the demand for it oc- 
curs, then that demand is assumed to be filled from 
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outside the RSS and the shortage penalty cost reflects 
the time required to obtain a unit from the external 
supply system. 

2. The procurement leadtime for orders of replenishment 
stock is assumed to be zero. This permits the optimal 
stock quantity computed during the review cycle at the 
end of one period to be immediately on hand at the be- 
ginning of the next period. This assumption is fairly 
realistic for items carried at the NARF which are also 
maintained in stock at the supporting NSC. 

3. The unit rates for ordering, surplus and shortage 
costs are treated as linear. Thus, for example, a 
cumulative excess of supply over demand of 3 units at 
the end of a particular period will result in a sur- 
plus penalty which is three times higher than if the 
total excess quantity was only 1. 



C . COST ELEMENTS 

There are many types of costs associated with maintaining 
inventory systems, ranging from obvious factors such as pro- 
curement and storage costs to less evident ones like the cost 
of operating data gathering procedures for the system [Ref. 

4: 10]. However, the only costs which need to be considered 

in calculating the optimal inventory quantity for a particu- 
lar item are those which vary with the quantity ordered. The 
components of the total variable cost (TVC) associated with 
the model developed in this thesis are discussed in the fol- 
lowing paragraphs. 

1 . Ordering Costs 

A basic assumption behind the MRP concept is that an 
order review will be conducted once per production period to 
determine if a stock replenishment action for each individual 
line item should be accomplished. Since this review will be 
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conducted at the close of each period for all items and since 



subsequent requisitions will be created as part of the review 
process, the cost of accomplishing the review and submitting 
the requisitions will be considered constant from period to 
period. Thus, the only variable cost associated with the or- 
dering process will be the cost of the quantity procured. 

This is the product of the unit price for a particular item, 
C, and the quantity of that item which is ordered. 

2 . Surplus Costs 

Surplus costs are those associated with having a 
quantity of an item on hand in inventory beyond the period in 
question. They can thus be expressed as the cumulative ex- 
cess of supply over demand at the conclusion of a given 
period times a constant surplus cost rate per item, H. Since 
the costs of storing RSS stock remain relatively unchanged 
from period to period, the surplus cost rate can be treated 
as a penalty paid whenever the supply system is unwilling to 
take back NARF excess stock on a full credit basis. This 
penalty can be quantified by combining the unit price of the 
item in question with the historical percentage of excess 
NARF material accepted by the supply system for credit. Ac- 
cordingly, surplus cost rates which are less than a particu- 
lar item's unit price will be utilized in evaluating specific 
model results. 
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3 . Shortage Costs 



Whenever a rework requirement for a spare part cannot 
be filled from the RSS inventory, a shortage cost represent- 
ing the penalty associated with a work stoppage is incurred. 
This shortage cost rate per item, P, is assessed at the con- 
clusion of a period on the cumulative excess of demand over 
supply. The shortage cost is typically the most difficult 
parameter in an inventory model to quantify since a reasona- 
ble estimate must include an accumulation of all production- 
related delay costs such as the labor charge to backrob or 
place the end item in storage until the required part arrives 
[Ref. 2: 32]. As a result, particular emphasis will be 
placed on examining the sensitivity of results obtained from 
specific models to changes in the shortage cost rate. 

D. SINGLE-PERIOD OBJECTIVE FUNCTION 

By utilizing the Karlin inventory model as adapted to in- 
clude particular NARF costs, an objective function expressing 
TVC for a single-period model can be formulated. This func- 
tion can then be used to identify that inventory order quan- 
tity for a particular item which results in the lowest TVC 
for a given demand probability distribution. As previously 
noted, only those items which are replaced less than 100% of 
the time during overhaul need to be examined in the model. 

Using the surplus and shortage cost parameters as pre- 
viously discussed, the expected costs at the end of one 
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period given an initial quantity of stock, y, and a continu- 
ous demand probability distribution, g(u), may be expressed 
as L (y ; g) , where : 

y 00 

L(y;g) = H J* (y-u)g(u)du + Pf (u-y)g(u)du. (1) 

0 y 

On the interval where demand, u, takes on values from 0 to y, 
the amount (y-u) represents a surplus quantity. Similarly, 
the amount (u-y) is a shortage quantity whenever u ranges 
from y to infinity. Therefore, the first term of L(y;g) is 
the product of the surplus cost rate per unit and the ex- 
pected number of surplus units at the end of the period, 
while the second term consists of the product of the shortage 
cost rate and the expected number of units short at the end 
of the period. 

The TVC incurred when x is the amount of an item on hand 
at the conclusion of the previous period and y is the quanti- 
ty in stock at the beginning of the present period can be 
calculated as: 

TVCy = C[y-x] + L(y;g) 

where C[y-x] represents the total procurement costs and 
L(y;g) was given by equation (1) . The minimum TVC for a one- 
period model will be denoted f(x;g), where: 

f(x;g) = min{C[y-x] + L(y;g)}. (2) 

ylx 



19 



The optimal y quantity used in identifying f(x;g) will be 
abbreviated as y*. 



E. DEMAND PROBABILITY DISTRIBUTIONS 

As discussed in the previous section, the value of y* 
obtained from a specific model is heavily dependent upon the 
applicable demand probability distribution. A realistic ex- 
pression of actual demand experienced at a NARF may be de- 
rived by using the binomial distribution. Under this 
approach, a particular repair part contained in a given* 
component will either be required or not required in any re- 
work action, and the probability, p, of the part being re- 
quired can be estimated using historical data available at 
the NARF. Given that the production schedule for a particu- 
lar component is n during the period in question, and that 
one unit of a certain repair part is contained in that com- 
ponent, the probability, g(u), that total demand for that 
repair part during this period will be u units may be ex- 
pressed as [Ref. 5: 84], 



g (u) = 



nl 



u! (n-u) 



u,T .n-u 
P (1-p) 



Recursion equations may be utilized to simplify the pro- 
cess of computing the entire range of demand for this dis- 
crete function. Thus, it can be easily shown that: 
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(3) 



g (u) 




[n-(u-l) ] 
u (1-p) 



(1-p) ^ for u = 0 



pg (u-1) for u > 0 



The binomial demand distribution as used above is appro- 
priate only for the situation in which a quantity of one of 
the repair part in question is present in the parent equip- 
ment being overhauled. The case involving more than one of 
a specific part in each of identical end items can also be 
shown to be binomial with an additional parameter to describe 
the multiplicity of the part [Ref. 2: 4], However, when a 

common item is present in several different equipments, its 
demand distribution is not binomial and will not be con- 
sidered further in this thesis. 

The binomial distribution is cumbersome to use in identi- 
fying optimal inventory quantities through the Karlin model. 
Thus, in order to illustrate the optimization process, the 
initial problem formulation will utilize the uniform distri- 
bution to represent predicted demand. Both continuous and 
discrete uniform distribution functions will be considered. 

Assume that the production schedule for a particular com- 
ponent during the period in question is n, and that one unit 
of a specified repair part is present in that component, but 
that not all of it needs to be replaced. An example of this 
would be adding additional oil to a partially full hydraulic 
reservoir. Then, all values of demand, u, between 0 and n 
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are equally likely under the continuous uniform distribution 
[Ref. 6: 113-114]. By restricting the range of possible de- 

mand values to integer quantities, the discrete uniform dis- 
tribution may be applied. In this situation, the probability, 
g(u), that total demand for the repair part in question will 
be u units, where u = 0 , 1 ,2 , . . . ,n, is given in general by: 

{ — for 0 i u ^ n 

(4) 

0 otherwise. 

F. SINGLE-PERIOD DERIVATIONS 

In this section, the specific single-period models asso- 
ciated with the continuous and discrete uniform demand dis- 
tributions as well as the binomial density function will be 
derived. For all models, the optimal stock quantity, y*, at 
the beginning of the period, given a production schedule of n 
units and a previous period ending balance of x, will be that 
value of y which minimizes the expected TVC . The model uti- 
lizing a continuous uniform demand distribution function will 
allow the y* value to be calculated explicitly using the cal- 
culus. However, the models based on discrete distributions 
require the approach of finite differences which involves 
comparison of TVC values associated with integer inventory 
quantities . 

The optimal inventory value computed by applying each 
model should be viewed both as a critical number which acts 
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as a cutoff point for determining whether an order should be 
placed as well as a high limit defining what the order quan- 
tity should be. Specifically, if y* is less than or equal to 
the quantity, x, on hand at the close of a given period when 
the level review is accomplished, then the optimal policy is 
not to order additional stock. Conversely, a y* value great- 
er than the closing inventory balance implies that a requisi- 
tion for (y*-x) is needed to bring available stock up to the 
identified starting optimal quantity for the period. 

1 . Continuous Uniform Distribution 

In this situation, let g(u) denote the continuous 
uniform demand distribution function whose characteristics 
have been previously described. Then, equation (2) may be 
directly used as follows in evaluating the optimal inventory 
quantity : 

f(x;g) = min{C[y-x] + L(y;g)} 
y^x 

y 00 

= min|c[y-x] + hJ* (y-u)g(u)du + P J' (u-y) g (u) du} . 
y-x ^ 

■^0 y 

Since the function g(u) is zero for values of u > n from equa- 
tion (4), the TVC function to be minimized may be rewritten 
as : 



y n 

TVC^ = C[y-x] f y* (y-u)du (u-y)du. 

0 y 
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Carrying out the integration results in: 



TVCy = + (C-P)y + - xC. 

Since Karlin has proven that the general TVC function re- 
sulting from his model is convex [Ref. 3: 236], minimization 

of this continuous function may be accomplished by taking the 
derivative with respect to y, setting it equal to 0, and 
solving for the optimal y. Thus, 

0 = + (c-P) 

dy n 

or. 



It can be seen from equation (5) that if the shortage penalty, 
P, is less than or equal to the unit cost, C, the optimal 
course of action is not to stock the item at all. This re- 
sult occurs in the P i C situation because it is cheaper to 
incur the penalty cost for being out of stock than it is to 
make the ordering cost investment needed to bring the item 
into the inventory. It should also be observed that the pre- 
vious period closing inventory balance, x, has no impact on 
the value of y* . 

2 . Discrete Uniform Distribution 

As previously noted, discrete demand distributions 
require a slightly different approach for evaluating y* than 
used with continuous functions because the required demand 
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probabilities are only defined for integer u values. Ac- 
cordingly, the method of taking finite differences must be 
applied. The technique of finite differences calls for de- 
termining the largest integer value of y such that the ex- 
pression (TVC - TVC .) is less than zero. In this manner, 
y y-1 

that initial stock quantity which generates the lowest ex- 
pected TVC will be identified since the cost function is con 
vex and it can be shown that non-integer y values can never 
be optimal. 

Assume that g(u) represents the discrete uniform de- 
mand distribution. To use this particular distribution in a 
single-period model, the TVC function must be modified to ac 
commodate the discrete case. This may be accomplished as 
follows; 

TVC^ = C[y-x] + L(y;g) 

y_l 00 

= C[y-x] + H ^ (y-u)g(u) + P I (u-y)g(u). 
u=0 u=y 

Substitution of g(u) from equation (4) results in: 

H P ^ 

TVC = C[y-x] + - I (y-u) + - I (u-y) . 

^ ^ u=0 ^ u=y 

By making the necessary substitutions and performing appro- 
priate cancellations, it can be shown that; 

TVC - TVC = y ^ . 

y y-1 n 
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Hence, y* is the largest y value such that the above formu- 
lation is negative. This condition corresponds to the lar- 
gest value of y for which: 

^ (n+l)P - nC 

y < ifpp • (6) 

As with the continuous case, the stockout penalty charge must 
be greater than the unit price if y* is to take on positive 
values. Also, y* is again independent of x. 

3 . Binomial Distribution 

The binomial distribution is a more complex example 
of a discrete demand function than its uniform counterpart 
since this distribution will take on different probability 
values for possible demand quantities within the production 
range. Again, the method of finite differences is appro- 
priate for determining y*. 

Accordingly, let g(u) represent the probability that 
demand during the period in question will be u units, given 
a production schedule, n, and a repair part replacement fac- 
tor, p. Then, as in the discrete uniform case, 

y-1 n 

TVC = C[y-x] + H ^ (y-u)g(u) + P 5! (u-y)g(u) 
y u=0 u=y 

y-1 y-1 n n 

= C[y-x] +Hy ^ g(u) - H ^ ug(u) + P ^ ug{u) - Py I g(u). 

u=0 u=0 u=y u=y 

By combining terms and recognizing that the mean of the bino- 
mial distribution is equal to np, this may be simplified to: 
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y-1 y-l 

TVC = C[y-x] + (H+P)y ^ g(u) - (H+P) ^ ug (u) + P (np-y) . 

^ u=0 u=0 



Using this result, it can be shown that: 



TVC„ - TVC„ T = (H+P) J g(u) + (C-P) . 
^ u=0 



As with the discrete uniform case, y* is the largest y value 
for which the previous expression is negative. This condi- 
tion corresponds to the largest value of y for which: 



i g(u) 

u=0 



P- 
H+P * 



(7) 



The same conclusions about the relationship between P and C 
as well as between y* and x can be drawn from this model as 
were identified in the two previous derivations. 



G. SINGLE-PERIOD EXAMPLE 

The following example is provided to illustrate the 
single-period models. Assume that the following data is 
available : 

n = 10 units, C = $100, H = $50, P = $200, x = 4 units. 



Using the continuous uniform density function to represent 
the demand distribution, the optimal opening inventory 
balance may be computed from equation (5). Thus, 



10 ( 200 - 100 ) 
50+200 



4 units. 
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since this quantity is equal to x, no order for additional 
stock should be placed. 

Similarly, if g(u) is assumed to be the discrete uniform 
demand distribution, equation (6) shows that y* will be the 
largest integer y value such that: 



Therefore, y* = 4 units. Again, the optimal procedure is not 
to order additional stock since a quantity of 4 units is al- 
ready available. 

Finally, in the binomial model with the replacement rate, 
p, assumed to be 0.5, equation (3) can be applied to generate 
the appropriate demand probability distribution values. This 
information can then be used in equation (7) to evaluate y* 
as 5 units. In this case, one additional unit should be or- 
dered to bring the current inventory balance up to y* . 

H. SUMMARY OF KEY FINDINGS 

The major results derived in this chapter as applicable 
to single-period inventory models may be summarized as 
follows: 

1. The optimization process identifies a critical number 
which is compared to the previous period closing 
balance in order to determine whether stock replen- 
ishment action is appropriate, and, if so, for what 
quantity . 

2. y* values may be obtained for the demand distributions 

examined by applying equations (5) , (6) , and (7) . 



y < 



( 11 ) ( 200 ) 




50 + 200 
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3. The shortage penalty charge, p, must exceed the unit 
price, C, for the models to produce optimal initial 
inventory values which are positive. 

4. y* is independent of the previous period closing 
balance, x. 
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III. DUAL-PERIOD MODELS — UNIFORM DISTRIBUTION 



A. DISCUSSION 

The single-period model developed in Chapter II will now 
be extended in an attempt to take advantage of production in- 
formation for an additional period. However, the requirement 
to project costs for the second period is dependent on the 
actual demand during the first period, thus making the two- 
period model considerably more complex than the single-period 
situation. The approach taken in this chapter will be to 
first describe the general form of the dual-period model as 
well as its multi-period counterpart. Then, the specific 
model associated with the continuous uniform demand distribu- 
tion function will be derived since this model permits expli- 
cit evaluation of the optimal starting inventory quantity. 
Finally, the effects of changing the demand distribution to 
the discrete uniform probability mass function will be ex- 
plored. This will establish a degree of familiarity in deal- 
ing with discrete demand functions and will lay the groundwork 
for the dual-period binomial model which is examined in Chap- 
ter IV. In all two-period system derivations, the general 
assumptions made in Section B of Chapter II will continue to 
apply . 
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B. OBJECTIVE FUNCTION — DUAL/MULTI-PERIOD MODELS 

As accomplished for the single-period case, the Karlin 
model can be adapted to provide an inventory system which de- 
termines optimal stock quantities based on projected produc- 
tion data for several quarters. This multi-period model 
assumes that demand is represented by a sequence of indepen- 
dent random variables which cover successive periods and are 
not necessarily identically distributed [Ref. 3: 233]. One 

familiar characteristic of such a model is that the resulting 
optimal inventory quantity for the first period is a critical 
number which has the same properties as did y* in the single- 
period system. An additional feature of the two-period model 
is that evaluation of the TVC function for any opening inven- 
tory value, y, is dependent on the computation of several 
optimal results from the second period considered alone. 

This will make determination of the optimal two-period inven- 
tory quantity a much more difficult task than was the case 
for the one-quarter model. 

1 . Dual-Period 

Assume that x and y represent the previous period 
closing inventory balance and the present quarter opening 
quantity, respectively. Additionally, let g and h be con- 
tinuous demand probability distributions for the next two 
periods. Then, the minimum expected TVC for a two-period 
model will be denoted f(x;g,h) where [Ref. 3: 235], 
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f(x;g,h) = min{C[y-x] + L(y;g) + f{0;h) J’ g (u) du 
y^x y 

Y 

+ (y-u;h) g (u) du} . (8) 

0 

As with the single-period case, that y value which generates 
f(x;g,h) will be called y*. 

In reviewing equation (8) , it should be noted that 
the first two terms correspond, respectively, to procurement 
and surplus/ shortage costs during the first period for the 
applicable opening inventory value. The third term repre- 
sents the expected optimal total costs incurred in the second 
period as a result of having no initial stock balance in that 
period. Finally, the last term is the expected total optimal 
costs for the second quarter, given a positive starting inven- 
tory balance in that period as a consequence of the proba- 
bility that demand during the first period is less than the 
initial stock quantity. 

2 . Multi-Period 

Although multi-period inventory systems beyond the 
dual-period case are not within the scope of this thesis, the 
general finite-period model as adapted from Karlin should 
also be mentioned [Ref. 3: 235]. Accordingly, assume that 

g^,g 2 ,...,g^ represent the demand probability distributions 
for the next n quarters. As with the preceding systems, the 
minimum TVC for the n-period model will be denoted 
f (x;g^,g 2 , . . . ,gj^) and evaluated as: 
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■m 











mm’ 




f (x;g^,g 2 , . . , ,gj^) =min{c[y-x] + L(y;g.) 



y-x 



CO 



+ f (0 ;g2 ,g3 , . . . ,g^) fq-^(u)dn 



y 



y 

+ y*f (y-u;g 2 ,g 3 , . . , ,g^^) gj^ (u) du} . 



0 



As can be seen, the individual terms in this formulation 
closely correspond to the previous description of a two- 
period model. 

C. GENERAL DUAL-PERIOD RESULTS 

Before proceeding to examine specific dual-period models, 
it is appropriate to summarize several important theorems 
generated by Karlin. The results which follow are dependent 
on the concept of stochastic ordering which is defined in the 
next paragraph. 

Assume that (g^ ,g 2 , . . . ,gj^) and (h^ ,h^ , . . . ,h^) represent 
two sequences of continuous demand distributions for periods 
1 through n. Then, the sequence g^ is said to be stochasti- 
cally smaller than the sequence h^ if: 



for i = l,2,...,n and all non-negative y values. In other 
words, demands based on the function g^ have a larger proba- 
bility of assuming smaller values than those generated by h^ . 
A simple single-period example of stochastic ordering can be 



y 



y 




0 



0 
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seen by comparing two continuous uniform distributions based 
on different production schedules. In this situation, the 
distribution resulting from the lower expected workload is 
stochastically smaller than the one associated with the high- 
er production schedule. 

Using this concept, several key results as proven by 
Karlin and then adapted to a dual-period model are listed as 
follows : 

1. If g and h are continuous demand distributions for the 
next two periods and g is stochastically smaller than 
h, then the starting optimal inventory quantity for the 
two-period model generated by the demand sequence (g,h) 
is greater than or equal to the single-period y* value 
determined solely by g. 

2. If g and h are continuous demand distributions for the 
next two periods and g is stochastically smaller than 
h, then the optimal inventory quantity for the first 
period based on the sequence (g,h) is equal to the y* 
value identified by the sequence (g,g). This is a par- 
ticularly useful result for verifying information pro- 
vided by specific models derived in the remainder of 
this thesis. 



D. CONTINUOUS UNIFORM MODEL 

The first specific dual-period model to be evaluated will 
be that system which is associated with the continuous uni- 
form demand probability distribution. Accordingly, let m and 
n be the production schedules for a particular component in 
two successive rework periods, and g and h be the correspond- 
ing continuous uniform demand distributions. By applying 
equation (8) , the optimal result for the two-period inventory 
model can be expressed as: 
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f(x;g,h) = min{C[y-x] + L(y:g) + f(0;h) rg(u)du 
y-x y 

y 

+ J' f (y-u;h) g (u) du} . 

0 

As will shortly be demonstrated, evaluation of the first 
three terms in this cost function is a relatively straight- 
forward process. However, two separate cases must be con- 
sidered in assessing the last term. 

1 . Common Objective Function Terms 

The individual terms in the cost equation may be 
simplified as follows: 

y m 

a. L(y;g) ~ (y-u)g(u)du + P J* (u-y)g(u)du. 

0 y 

Applying the appropriate value of g and carrying out the in- 
tegration results in: 



L (y;g) 



(H+P)y^ _ 
2m 



Py + 



mP 

2 



00 m 

b. J'g(u)du = J*g(u)du = . 

y y 

c. The previous chapter demonstrated the existence of a cri- 
tical number associated with any single-period continuous 
uniform model (see equation (5) ) . This number determines 
whether or not a stock reorder should be submitted as well as 
the appropriate quantity in those situations requiring re- 
order action. Accordingly, let k denote the critical number 
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pertinent to just the second quarter of this two-period 
model. Based on previous results, a first period ending 
balance of zero implies that a requisition for k units must 
be submitted to obtain the optimal inventory balance at the 
start of the second period. Thus, 

f(0;h) = kC + L(k;h) 

k 00 

= kC + h/ (k-u)h(u)du + p/ (u-k)h(u)du. 

0 k 

When this is evaluated, it reduces to: 

f(0;h) = + (C-P)k + ^ . 

d. Now, in the general case, the properties of k may be ap- 
plied to show that: 

rL(y-u;h) if y-u 1 k 
f(y-u;h) = < 

l^C[k-(y-u)] + L(k;h) if y-u < k. 

Therefore, it may be concluded that: 

y y-k 

J* f (y-u;h) g (u) du = f L (y-u ; h) g (u) du 
C 0 

y 

+ J" {c[k-(y-u)] + L (k;h) }g (u) du. (9) 
y-k 

In the case where y < k, the first term drops out while the 
lower bound over which the integration is performed on the 
second term becomes zero. Thus, two different cases need to 
be evaluated. 
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Case 1 



2 . 

In this case, y i k, so: 
y-k y-k 

J* L (y-u;h) g (u) du (y-u) ^ - P (y-u) + 

0 0 

Also , 

y T y 

J* {c tk- (y-u) ] + L (k;h) }g (u) du = — J' C[k-(y-u)]du 
y-k y-k 

y-k 

After evaluating these integrals, the entire objective func- 
tion for case 1 may be constructed by combining this result 
with the formulas for the common objective function terms 
from subsection 1 above. Finally, after taking the first 
derivative with respect to y and setting it equal to 0 , the 
following result is obtained: 

* _ -nH + Vn^H^-n[P^ (n-2m) +nC^+2CP (m-n) +2mH (C-P) ] 
y H + P • ^ ^ 

3 . Case 2 

Since y < k in this case, the last term in the objec- 
tive function given by equation (9) is now: 

y y 

J" f (y-u;h) g (u) du = J' {C [k- (y-u) ] + L (k ,h) }g (u) du 
0 0 

= If C(k-(y-u)ldu + iJ {i«^^-plc + np}du. 
0 0 
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After carrying out the integration, the entire cost function 
for case 2 is formed as before by adding this result to the 
common objective function terms. Using the calculus, the 
optimal y value can then be identified as; 



* ^ m(P-C) 

H+P-C • 



( 11 ) 



It should be noted that, in both cases, the formula for com- 
puting y* is independent of the closing stock balance, x, 
from the previous period. This result corresponds to an 
identical finding for the single-period continuous uniform 
model. Also, the value of P must again be greater than C for 
positive y* values to occur. 

4 . Breakpoint Between Cases 1 and 2 

Prior to examining the results obtained from this 
model for specific parameter values, it is useful to note 
that the two distinct cases are equivalent at y = k. There- 
fore, the breakpoint between these alternatives may be calcu- 
lated in terms of parameter values by setting the optimal 
result from either case equal to the expression for k as de- 
rived from equation (5) . In particular, using the case 2 
result , 



m(P-C) ^ n(P-C) ^ , 
H+P-C H + P 



Therefore , 



P 



nH - nC - mH 
m - n 



( 12 ) 
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For P values greater than the right-hand side of this ex- 
pression, equation (11) should be used to calculate the op- 
timal inventory quantity, while smaller P values require the 
use of equation (10). Equation (12) may also be utilized to 
obtain a definition of this breakpoint between the cases in 
terms of the other parameters. 

E. EXAMPLE OF THE CONTINUOUS UNIFORM MODEL 

Assume that the projected production schedule for the 
next two quarters is 100 and 200 units, respectively. Addi- 
tionally, let C = $30, H = $20 and x = 0. Using the previous 
findings, a comparison of the optimal stock quantities at the 
beginning of the next period for the single and dual-period 
models may be conducted for various values of P. The single- 
period optimization utilizes the projected first period work- 
load but ignores subsequent production information. Curves 
of the y* values obtained from both models are displayed in 
Figure 3.1. In the dual-period model, a P value of 40 was 
computed as the breakpoint between cases 1 and 2 by using 
equation (12) . 

These results are intuitively appealing for several rea- 
sons. First, increasing the stockout cost causes the optimal 
initial inventory to take on larger and larger values under 
both systems. Additionally, the dual-period model, which 
considers projected demand beyond the first period, computes 
constantly larger inventory quantities than the single-period 
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Figure 3 . 1 



Continuous Uniform Model: One vs Two-Period 



model. Since the demand distributions are stochastically in- 
creasing, this is in agreement with the findings of Karlin 
which were previously discussed. Finally, extremely large 
stockout penalty costs generate nearly equivalent results 
under both systems as the initial inventory is maintained at 
a high level to protect against this exorbitant penalty 
charge . 

An assumption in this example was that the production 
schedule in the second quarter exceeded that of the first 
period. Suppose that, instead, the rework requirements are 
projected to be decreasing over time. Let the parameters 
used in the first example remain the same, except that now 
m = 200 and n = 100. The curves of y* corresponding to the 
single and dual-period models are shown in Figure 3.2. The 
y* values for the dual-period system are derived entirely 
from case 2 since the application of the breakpoint formula 
yields a negative value of P. 

These values are similar to the results of the increas- 
ing production schedule situation as the two-period system 
generates higher optimal inventory quantities than the one- 
quarter model despite the projected decrease in rework ac- 
tivity. This leads to the general conclusion that a 
two-period formulation can be applied in both increasing and 
decreasing production schedule situations to identify a y* 
solution yielding lower TVC results than those derived from a 
single-period y* value. 
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Two-Period 




Figure 3.2 



Continuous Uniform Model: Decreasing Production 
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F. DISCRETE Ul'IIFORJI MODEL 

A logical extension, of the continuous uniform model is 
to consider the discrete case in which both demand and in- 
ventory quantities may take on only integer values. This 
formulation provides a closer representation of the actual 
situation at the NARF where integer values of these parame- 
ters are appropriate. Since the demand function under this 
model is not continuous, finite differences are needed in 
determining y*. 

As before, assume that the production schedules for the 
next two periods are m and n respectively. Equation (8) can 
then be modified to accommodate' the discrete case as follows: 



r ra 

f(x;g,h) = minjc[y-x] + L(y;g) + f(0;h) I g(u) 
y^x I u=y 

y-1 1 

+ ^ f (y-u;h) g (u) k 

u=0 J 

Now, TVCy - TVC^_^ may be examined on a term by term basis as 
follows: 

1. C[y-x] - C[ (y-1) - x] = C. 

y-1 m 

2. L(y;g) - L(y-l;g) = (y-u)g(u) + (u-y)g(u) 

u=0 u=y 

y-2 m 

- [ (y-1) -u]g(u) - P [u- (y-1) ] g (u) . 

u=0 u=y-l 

After several steps, this may be reduced to: 

. , . V (H+P)y - (m+l)P 

L(y;g) - L(y-l;g) = . 
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3. As before. 



f(0;h) = kC + L(k;h) 

H ^ 7 ^ P ^ 

=kC + - I (k-u) + ^ I (u-k) 

" u=0 " u=k 

which is independent of y. The finite differences result for 
the third term is therefore 




4. The last term of the cost equation must again be con- 
sidered in two separate cases. For case 1, 

y-1 y-k y-1 

^ f (y-u;h) g (u) = ^ L (y-u; h) g (u) + C [k- (y-u) ] g (u) 

u=0 u=0 u=y-k+l 

y-1 

+ Z. L (k;h) g (u) . 
u=y-k+l 

However, the application of the method of finite differences 
in this case fails to provide the necessary cancellations 
which simplify the evaluation process and, thus, justify use 
of this technique. As a result, a complete enxameration of 
the TVC results for all possible y values offers a better ap- 
proach to solving case 1 than does finite differences. Since 
this procedure is employed in Chapter IV to evaluate the 
dual-period binomial model, it will not be duplicated here. 
Rather, the results for case 2 under finite differences will 
be shown . 



44 



Case 2 



Since y < k. 



y-1 y-1 y -1 

I f(y-uih)g(u) = I C [k-> (y-u) ] g (u) + I L(k;h)g{u). 
u=0 u=0 u=0 



Taking finite differences results in: 



= C(k2i ^ ^ X I 

m mn mn 



u=0 



u=k 



Putting the previous results all together, the following 
expression for the difference in variable costs between 
quantities of stock y and (y-1) may be obtained for case 2 
of the two-period discrete uniform model: 

mvr _ Tvr = (H+P-C)y + mC - (m-H)P 

ivuy m 



Thus, the optimal y is the largest integer quantity for which 
the above expression is negative. This condition corresponds 
to the largest value of y for which: 



y < 



(m+1) P - mC 
H+P-C 



(13) 



If there is also an integer value of y which is equal to the 
right-hand side of equation (13) , then it is an alternate 
optimal solution. As in the previous model, the optimal y 
value is independent of the previous period closing balance, 
X, and P must be greater than C for y* to be positive. 



45 



G. EXAMPLE OF THE DISCRETE UNIFORM MODEL 

Figure 3.3 presents the results obtained for a range of 
stockout costs utilizing the discrete uniform demand distri- 
bution in single and dual-period models as derived above and 
in Chapter II. The values for m, n, x, C and H remain the 
same as in Figure 3.1. Only relatively large values of P 
are analyzed to ensure that the breakpoint between cases 1 
and 2 is exceeded. The y* values obtained are consistent 
with those computed from the continuous uniform approach. 



H. SUMMARY OF KEY FINDINGS 

The principal results obtained in this chapter for two- 

w 

period inventory models based on uniform demand distributions 
are sxammarized as follows: 

1. As in the single-period case, a critical number which 

determines the necessity of ordering additional stock 
at the beginning of the first period as well as the 
optimal order quantity can be identified. This number 
is independent of the previous period closing balance, 
and may be computed for the various cases described by 
applying equations (10) , (11) and (13) . Additionally, 

the shortage cost rate must be greater than the unit 
price if the item is to be stocked at all. 

2. The critical number for the first quarter of a dual- 
period model is always greater than or equal to the 
optimal result for the corresponding one-period model. 

3. In both increasing and decreasing production schedule 
situations, the two-period model identifies a y* quan- 
tity which yields lower total costs than those gen- 
erated by a one-period y* value. 
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Figure 3.3 Discrete Uniform Model: One vs Two-Period 
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IV. DUAL-PERIOD BINOMIAL MODEL 



A. DISCUSSION 

Now that an understanding of the general rationale behind 
a two-period inventory system has been developed using the 
uniform demand distribution, attention will turn to the bi- 
nomial density function which is more pertinent to a NARF. 

As with the discrete uniform case, explicit computation of 
the optimal inventory quantity for this dual-period model 
will not be possible since the demand distribution is only 
defined at integer values. Therefore, a technique of com- 
paring the TVC values for possible initial inventory quanti- 
ties must be employed. 

The computation of total costs for this distribution will 
be more complex than for the uniform distribution since the 
density function will take on different values for various 
demand possibilities within the production range. The same 
complication will also be encountered as in the previous ex- 
amples whereby separate evaluation of the last term of the 
objective function for two distinct cases must be accom- 
plished. As a consequence, the technique of finite differ- 
ences will not provide any advantages over using the TVC 
values directly. Accordingly, the general approach to this 
particular model will be to first derive the entire TVC for- 
mulation for both cases, and then to use these formulas to 
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evaluate the expected variable costs incurred for each possi- 
ble starting inventory value. From this enumeration, the 
optimal y associated with the minimum TVC can be selected. 

B . THE MODEL 

Let g and h represent the binomial demand distributions 
for the next two periods derived from an historical replace- 
ment factor, p, and production schedules m and n. As with 
the discrete uniform model, equation (8) must be modified to 
yield an expression for TVC associated with the opening in- 
ventory balance, y, and a previous period ending quantity, x. 
The y* value has been shown to be independent of the previous 
period closing balance for all models previously evaluated in 
this thesis, and a similar rationale can be employed in this 
particular situation. Therefore, x will be assumed to be 
zero for convenience in formulating the model. Then, 

m y-1 

TVC = Cy + L(y;g) + f(0;h) ^ g(u) + I f (y-u; h) g (u) . 

^ u=y u=0 

1 . Common Objective Function Terms 

The evaluation of the first three terms of the objec- 
tive function will be accomplished independently of the last 
term which involves separate cases contingent on the value of 
the critical number for the second period. Thus, 
a. The first term, Cy, and one component of the third term, 
g(u), cannot be simplified. 
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y-1 m 

b. L(y;g) = H ^ (y-u)g(u) + P ^ (u-y)g(u) 
u=0 u=y 

y-1 y-l m m 

= Hy ^ g(u) - H J ug(u) + P ^ ug(u) - Py 5! g(u) . 
u=0 u=0 u=y u=y 

However, the following substitution may be made: 

m m y-1 y-1 

I g(u) = 1 g(u) - \ g(u) = 1 - I g(u) . 

u=y u=0 u=0 u=0 

Also , 

m m y-1 

I ug(u) = I ug(n) - I ug(u). 
u=y u=0 u=0 

Now, the first term on the right-hand side of the expression 
immediately above is the mean of g(u). Hence, 

m y-1 

I ug (u) = pm - 5! • 

u=y u=0 

With these substitutions, L(y;g) reduces to: 

L(y;g) = (H+P)y I g(u) - (H+P) I ug (u) +pmP-Py. (14) 

u=0 u=0 



c. As in previous formulations, let k denote the critical 
number derived from an independent evaluation of the second 
period. This may be accomplished by using equation (7) as 
developed in Chapter II. Then, using the variable s to repre- 
sent demand in the second period, 

f(0;h) = kC + L(k;h) 

k-1 n 

= kC + H ^ (k-s)h(s) + P ^ (s-k)h(s). 
s=0 s=k 



50 



Applying the same simplifications as were used to obtain 
equation (14) results in; 



k-1 k-1 

f(0;h) = (C-P)k + (H+P)k i h(s) - (H+P) I sh(s) + pnP. 

s=0 s=0 

2 . Case 1 

For the case where y ^ k, we recall from Chapter III 

that: 

y-1 y-k y-1 

^ f (y-u;h) g (u) = L (y-u;h) g (u) + ^ C [k- (y-u) ] g (u) 

u=0 u-0 u=y-k+l 

+ ^ L (k;h) g (u) . 

u=y-k+l 



Now, this expression can be simplified using the definition 
of L(y;h) as contained in equation (1), the properties of the 
binomial distribution, and the variable s to represent demand 
in the second period. After simplification, the entire ob- 
jective function for this case becomes: 



y-^ 

TVC = (C-P) y + (H+P-C)y X, - (H+P-C) X, 

y u=o u=o 

y-k y-k 

+ pmP + Cy g (u) - C V ug (u) 
u=0 u=0 

m k-1 k-1 

+ Yi j(H+P)k ^ h(s) - (H+P) ^ sh(s) + pnP + k(C-P)[ g(u) 
u=y-k+l s=0 s=0 

y-k f y-u-1 y-u-1 

+ ^ j (H+P) (y-u) Y h(s)-(H+P) Y sh (s) + pnP - P (y-u)[g (u) . 
u=0 ^ s=0 s=0 

(15) 
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3 . Case 2 



In this situation with y < k, 

y-1 

I f(y-u;h)g(u) = I C [k- (y-u) ] g (u) + I L(k;h)g(u). 
u=0 u=0 u=0 

Using an identical approach to that of case 1, an evaluation 
of this expression may be performed. After retrieving all 
terms in the objective function, making the appropriate com- 
binations, and noting that the summation from 0 to m of g(u) 
equals 1, the result for case 2 is: 

y-1 y-1 

TVC = (C-P)y + (H+P-C)y ^ g(u) - (H+P-C) I ug(u) + pmP 
^ u=0 u=0 

k-1 k-1 

+ (H+P)k i h(s) - (H+P) I sh(s) + pnP + (C-P)k. (16) 

s=0 s=0 

4 . Conclusion 

The critical number y* for a two-period binomial 
model may thus be identified by using equations (15) and (16) 
to compute the TVC values for all possible initial inventory 
quantities; it is that value of y which results in the mini- 
mum TVC. It should be noted that, if the assumption of x 
equaling zero is dropped, then the TVC equation for each y 
value will include a Cx term. Therefore, the relationship 
between TVC values associated with all y quantities will be 
unchanged by the addition of this common term. 
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C. EXAMPLES OF THE BINOMIAL MODEL 



Because the objective functions developed in the two- 
period binomial model are rather complex, a computer program 
was written which accepts various system parameters and de- 
termines the optimal y value for the beginning of the first 
period. Unlike the continuous uniform model, identification 
of the breakpoint between cases 1 and 2 is not required as 
the program first determines k and then uses this value in 
calculating appropriate variable costs. 

Results obtained from using this program for various 
values of P with 

m = 10, n = 20, p = 0.5, x = 0, H = $20, C = $30 

are displayed in Figure 4.1. As in all previous single and 
dual-period models, values of P less than C imply that the 
item should not be stocked at all. Consequently, only P 
values above $30 are included in this figure. 

These results are consistent with those obtained using 
the uniform distribution models in that the optimal value of 
y at the beginning of the first period increases with the 
stockout penalty. Additionally, the extra information con- 
tained in the two-period model causes its optimal result to 
be greater than or equal to that obtained from the single- 
period analysis. However, the difference in the y* values is 
never more than one. Thus, an examination of the benefits 
derived from the additional stock quantities in the dual- 
period model needs to be conducted. 
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Figure 4.1 Binomial Model: One vs Two-Period 
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However, before performing this detailed analysis, it is 
useful to examine a decreasing production schedule situation 
as well as verify Karlin's result that a two-period inventory 
model based on stochastically increasing demand distributions 
generates y* values which are equal to the y* quantities pro- 
duced by a dual-period inventory system based on a constant 
production schedule. Accordingly, assume that the following 
parameter values are given: 

m = 10, X = 0, H = $20, C .= $30, P = $150 

The optimal initial inventory quantities derived for various 
replacement rate values and second-period production sched- 
ules are shown in Table I. 



TABLE I 



Binomial Model: Varying Production 



y* (n = 5) 



y* (n=10) 



(n = 20 ) y* (n ^ 50) 



0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 



10 

9 

8 

8 

6 

5 

4 

3 

2 



10 

9 

9 

8 

7 

6 

5 

3 

2 



10 

9 

9 

8 

7 

6 

5 

3 

2 



10 

9 

9 

8 

7 

6 

5 

3 

2 



Several conclusions 
the situation where the 



may be drawn from Table I. First, 
production schedule in the second 
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period is greater than or equal to ten fulfills the condi- 
tions for a stochastic demand increase, and, as noted by 
Karlin, this environment always produces a y* value equal to 
that of the constant production schedule situation. Addi- 
tionally, the case where demand is stochastically decreasing, 
as depicted in the table under n = 5, identifies y* values 
which are nearly equivalent to the increasing production 
schedule situation. This may be interpreted to mean that the 
mere prospect of any demand in the second period is more im- 
portant than the actual forecast quantity. 



D. SUI-U-IARY OF KEY FINDINGS 

The major results of this chapter for the two-period bi- 
nomial inventory model are summarized as follows: 

1. The optimal initial inventory value is again a critical 
number which may be computed by identifying that y 
value which generates the minimum TVC in equations (15) 
and (16) . As in the previous single and dual-period 
models, y* is independent of the previous period closing 
balance, x, and will only take on positive values if 
the shortage cost rate for an item exceeds its unit 
price . 

2. An increasing production schedule will produce a y* 
value identical to the situation in which the produc- 
tion workload is constant. 



56 



V. TWO-PERIOD BINOMIAL MODEL SENSITIVITY ANALYSIS 



A. DISCUSSION 

In this chapter, a detailed analysis of the two-period 
binomial model is conducted. The first portion of this 
analysis is a parametric evaluation of the sensitivity of 
model results to changes in the unit price and the surplus 
and shortage cost rates. Then, a comparison of expected 
optimal TVC values generated by several alternative binomial 
models using a two-period time horizon is conducted. These 
alternative costs are arrived at by employing the standard 
two-period model, the two-period objective function using 
single-period y* values, and the single-period model for two 
distinct periods. Based on the results of this analysis, 
several general conclusions on the use of the two-period bi- 
nomial model are drawn. 

B. ORDERING, SURPLUS AND SHORTAGE COST SENSITIVITY ANALYSIS 

One aspect of analyzing the dual-period binomial inven- 
tory model is to evaluate the sensitivity of the model to 
changes in various system parameters. Since the production 
schedules which generate the demand distributions are deter- 
mined outside the scope of the model, the current emphasis 
will be on examining the effects of changes in the unit price 
and surplus/ shortage cost parameters on system results. 
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Throughout this section, a constant production schedule of 10 
units for each period is assumed, while the replacement fac- 
tor is permitted to take on the values of 0.9, 0.5 and 0.1. 
Then, two of the three parameters of interest will be held 
constant so that the results obtained from varying the third 
may be examined. 

1 . Unit Price Changes 

An assumption was made in Chapter II that the surplus 
cost rate for an item must be less than its unit price. Ad- 
ditionally, the preliminary analysis conducted on the bino- 
mial model in Chapter IV established that the shortage cost 
rate should be greater than the unit price for any particular 
part to be stocked. Thus, a range of potential values for C 
between H and P seems appropriate. Under the assumption that 
H = $250 and P = $1000, Table II displays optimal initial in- 
ventory quantities associated with the two-period model and 
the previously identified replacement factors for various 
values of C within the established range. The TVC calculated 
at y* and at values of y which differ from y* by one unit are 
also shown. 

The general conclusion that may be drawn from this 
table is that y* is not particularly sensitive to changes in 
the unit price, except when the value of C approaches the 
ceiling level established by the shortage cost rate. How- 
ever, when C does get close to P, the difference between 
expected variable costs for y* and the initial inventory 
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TABLE II 



Binomial Model; Sensitivity Analysis of C 



(1) p = 0.9 



c 


y * 


TVC (y*-l) 


TVC (y*) 


TVC (y*+l) 


250 


10 


5284.52 


5185.84 





500 


10 


9697.35 


9685.86 


— 


750 


9 


14148.34 


14030.29 


14105.94 


950 


8 


17322.59 


17293.65 


17322.81 


990 


7 


17882.69 


17876.01 


17884.26 


p = 0.5 


C 


Y * 


TVC (y*-l) 


TVC (y*) 


TVC (y*+l) 


250 


6 


3884.28 


3757.35 


3836.04 


500 


6 


6230.46 


6198.02 


6322.36 


750 


5 


8416.99 


8355.68 


8419.88 


950 


3 


9818.16 


9784.58 


9786.48 


990 


2 


9981.47 


9974.27 


9978.57 


p = 0 . 1 


C 


y * 


TVC (y*-l) 


TVC (y*) 


TVC (y*+l) 


250 


1 


2000.00 


1284.53 


1329.94 


500 


1 


2000.00 


1697.36 


1895.92 


750 


1 


2000.00 


1990.08 


2413.22 


950 


0 


— 


2000.00 


2189.28 


990 


0 


— 


2000.00 


2230.08 



quantities immediately above and below y* is so small as to 
be insignificant. In fact, if a y value of 6 is used in the 
p = 0.5 model with C = $990, the expected TVC is only three 
percent higher than the corresponding TVC associated with the 
optimal y value of 2. 

2 . Surplus Cost Rate Changes 

In this portion of the analysis, the range of poten- 
tial surplus cost values will not be assumed to be bounded by 
other parameters, despite the assumption in Chapter II that 



59 



the surplus cost rate would be expected to be less than the 
unit price. This assumption about H is temporarily suspended 
so that the results obtained from a wide range of H values 
may be examined. Dual-period binomial model results for 



various H values were obtained for the three replacement fac- 
tors when C = $500 and P = $1000, and these values are listed 
in Table III. 

TABLE III 

Binomial Model: Sensitivity Analysis of H 

(1) p = 0.9 



H 


y* 


TVC (y*-l) 


TVC (y*) 


TVC (y*+l) 


25 


10 


9540.45 


9382.39 


— 


100 


10 


9592.75 


9483.54 


— 


250 


10 


9697.35 


9685.86 


— 


500 


9 


10107.78 


9871.69 


10023.00 


750 


9 


10216.14 


10046.01 


10360 . 17 


1000 


9 


10296.70 


10192.55 


10669.53 


1500 


8 


10741.10 


10381.48 


10409.28 


II 


0.5 

y* 


TVC (y*-l) 


TVC (y*) 


TVC (y*+l) 


25 


7 


5780.84 


5717.30 


5723.55 


100 


7 


5919.90 


5918.99 


5998.86 


250 


6 


6230.46 


6198.02 


6322.36 


500 


5 


6595.70 


6472.72 


6596.83 


750 


5 


6714.84 


6636.19 


6967.13 


1000 


4 


7076.17 


6833.98 


6899.66 


1500 


4 


7228.52 


7072.26 


7326.60 


II 


0.1 

y* 


TVC (y*-l) 


TVC (y*) 


TVC (y*+l) 


25 


2 


1540.45 


1515.64 


1760.36 


100 


1 


2000.00 


1592.75 


1642.40 


250 


1 


2000.00 


1697.36 


1895.92 


500 


1 


2000.00 


1856.89 


2312.86 


750 


1 


2000.00 


1974.51 


2712.50 


1000 


0 


— 


2000.00 


2092 . 12 


1500 


0 


— 


2000.00 


2327.35 
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The optimal initial inventory quantity decreases 
with increasing H but appears to be fairly insensitive to 
changes in the surplus cost rate. Additionally, the expected 
TVC at y* for H = $1500 are less than twenty-five percent 
higher than the results for H = $25 in all three cases, de- 
spite the fact that the surplus cost rate has been increased 
by a factor of sixty. These results imply that a precise 
estimate of the surplus cost rate is not required to obtain 
useful model results. 

3 . Shortage Cost Rate Changes 

Finally, the last parameter to be considered in this 
section of the analysis is the shortage cost parameter. As 
was noted in an earlier chapter, this parameter must be 
greater than the applicable item's unit price if the part in 
question is to be stocked at all. Thus, the range of values 
for P will start at C and has no upper bound. Table IV de- 
picts the dual-period model results obtained for various 
values of P when C = $500, H = $250 and the replacement fac- 
tor takes on the three values as before. 

The evidence from this table is that the optimal 
stock quantity at the beginning of the first period is in- 
sensitive to changes in the stockout cost rate when P is much 
larger than C. However, as the difference between the unit 
price of an item and its stockout cost grows small, the y* 
value changes more rapidly for relatively minor changes in P, 
especially in the case where p = 0.5. Additionally, the 
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TABLE IV 



Binomial Model: Sensitivity Analysis of P 



(1) p = 0.9 



p 


y* 


TVC (y*-l) 




TVC (y*) 


TVC (y*+l) 


550 


8 


9216.04 




9187.09 


9216.25 


750 


9 


9627.15 




9509.10 


9584.75 


1000 


10 


9697.35 




9685.84 


— 


2000 


10 


10360 . 17 




9999.98 


— - — 


5000 


10 


11406.20 




9999.98 


— 


10000 


10 


13149.43 




9999.99 


— 


20000 


10 


16636.25 




9999.99 


— 


(2) p = 0.5 












P 


y* 


TVC (y*-l) 




TVC (y*) 


TVC (y*+l) 


550 


3 


5306.64 




5273.04 


5274.69 


750 


5 


5857.42 




5796.02 


5859.15 


1000 


6 


6230.46 




6198.02 


6322.36 


2000 


7 


6953.12 




6902.70 


7061.15 


5000 


8 


7664.06 




7654.75 


7858.92 


10000 


8 


8328.12 




8045.11 


8192.78 


20000 


8 


9298.81 




8468.81 


8507.14 


(3) p = 0.1 












P 


y* 


TVC (y*-l) 




TVC (y*) 


TVC (y*+l) 


550 


0 


— 




1100.00 


1245.67 


750 


1 


1500.00 




1458.86 


1784.82 


1000 


1 


2000.00 




1697.36 


1895.92 


2000 


2 


2394.71 




2237.36 


2631.49 


5000 


2 


3851.30 




2847.77 


2897.38 


10000 


3 


3695.53 




3271.82 


3624.67 


20000 


3 


4719.71 




3583.43 


3751.84 


difference in 


expected TVC between 


y* and 


initial inventory 


levels immediately 


above and below 


y* becomes much greater as 


P takes on larger values. For 


example, the optimal result of 


y = 8 for P = 


$20000 and p = 0. 


5 is nearly 


ten percent less 


expensive than the 


expected TVC 


: at 


y = 7. 


The conclusion 


which may be 


drawn 


is that the 


NARF needs ■ 


to have reasonably 
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accurate estimates of its shortage cost rates, particularly 
for high value items. 



C. TVC COMPARISON — SINGLE VS DUAL-PERIOD MODELS 

In this section, the difference in the expected TVC be- 
tween single and dual-period binomial models is evaluated and 
analyzed. Comparison between the results obtained from the 
following distinct models is accomplished over a two-period 
cost horizon: 

1. The first alternative involves^ the computation of ex- 
pected total costs associated with y* as derived from 
the dual-period binomial model. 

2 . The second alternative computes an expected TVC value 
by a combined use of single and dual-period models. 

This is accomplished by substituting the y* value de- 
termined from a single-period model into the two-period 
objective function. This alternative, when contrasted 
with the pure dual-period result, depicts the differ- 
ences between one and two-period models on a scale 
which allows for a meaningful comparison of the systems. 

3. The third TVC result is determined by considering each 
quarter as a separate, single-period model. The y* 
values obtained for each period are then used in the 
single-period binomial objective function to calculate 
the aggregate expected TVC result. This alternative 
implies that any surplus at the end of the first quar- 
ter is disposed of and a buy for a quantity equivalent 
to the second-period y* is initiated. 

Assume that the model parameters take on the following 

values ; 

C = $500, H = $250, P = $1000, p = 0.5, x = 0. 

Table V displays the TVC results for the three models. The 
numbers in parentheses immediately following the total costs 
for each alternative identify the y* value (s) used in that 
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TABLE V 



Binomial Model: TVC Comparisons 



Production 

Schedules 

m=5, n=5 
m=5 , n=10 
m=5, n=20 
m=10, n=5 
m=10, n=10 
m=10, n=20 



Alternative 
One — TVC 

3316.41 (3) 
4808.11 (3) 

7625.41 (3) 
4747.93 (5) 
6198.02 (6) 
9015.06 (6) 



Alternative 
Two — TVC 

3437.50 (2) 
4933.11 (2) 

7750.41 (2) 
4747.93 (5) 
6230.46 (5) 
9047.78 (5) 



Alternative 
Three — TVC 

3546.87 (2,2) 
5042.48 (2,5) 
7859.79 (2,9) 
5042.48 (5,2) 
6538.08 (5,5) 
9355.39 (5,9) 



model in determining the expected TVC. The second of these 
numbers in parentheses for the third alternative is the k 
value for the second period which was described earlier in 
this thesis. 

As anticipated, the optimal inventory value identified 
by the two-period binomial model produces a lower expected 
TVC than either of the other two alternatives. The more 
expensive result for the second alternative results from the 
optimal dual-period y* not being used, while the increased 
costs associated with the third alternative are a conse- 
quence of disposing of the first period excess and, thus, 
having no carryover inventory. Despite these differences, 
the cost variations between the various model results are 
relatively small. Table VI depicts the differences between 
alternatives one and two as well as between one and three as 
a percentage of the expected TVC calculated for alternative 
one . 
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TABLE VI 



Binomial Model: Percentage Differences 



Production 

Schedules 



Alternatives 
One and Two 



Alternatives 
One and Three 



m=5, n=5 



3.7 

2.6 

1.6 



6.9 



m=5, n=10 
m=5, n=20 



4.9 



3.1 



m=10, n=5 



0.0 



6.2 



m=10, n=10 
m=10, n=20 



0.5 



0.4 



5.5 

3.8 



The most interesting result that can be seen from Table 
VI is that the percentage difference between alternative one 
and both other alternatives decreases as the production 
schedule for the second period increases, except in the case 
where the y* value for alternatives one and two is the same. 
In fact, by the time that the rework forecast for the first 
period reaches ten, the expected cost benefit to be derived 
from using the two-period y* instead of the single-period y* 
borders on the insignificant, thereby making the use of the 
more complex two-period model a questionable proposition. 

One additional factor which needs to be considered before 
concluding this model comparison is the accuracy of produc- 
tion schedule information. The first quarter workload fore- 
cast is utilized in both single and dual-period models, and 
so an error in this parameter will have a similar impact on 
the results obtained from each system. On the other hand, 
since the second period projected rework schedule is only 
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considered in the two-period model, a variance in this pa- 
rameter between the actual and forecasted workload will only 

m 

affect the results obtained from the two-period model. Addi- 
tionally, some difference between the actual and projected 
schedules is more likely to occur in the second period than 
in the first since current information is generally more ac- 
curate than that applying to a future period. While the im- 
pact of an inaccurate second period rework forecast will vary 
greatly from situation to situation, it will only serve to 
further reduce or even eliminate the limited advantages which 
the two-period model offers over its single-period counterpart. 

D. SUMMARY OF KEY FINDINGS 

The principal results identified in this chapter as ap- 
plicable to the dual-period binomial model are listed as 
follows : 

1. The model is very sensitive to the value of the short- 
age cost parameter and is relatively insensitive to the 
value of the surplus cost parameter. 

2. The cost advantages obtained from the two-period model 
represent only a small improvement over the one-period 
system. 

3. The difficulties in using the two-period model caused 
by its complexity as well as in obtaining accurate 
forecast data for the second quarter reduce that 
model's effectiveness. 
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VI. CONCLUSIONS /RECOMMENDATIONS 



A. DISCUSSION 

The basic purpose of this thesis has been to consider 
projected production information spanning two periods in com- 
puting RSS inventory quantities. One important assumption 
used throughout this thesis is that an order review is con- 
ducted at the beginning of each period for all RSS items, and 
that, consequently, there is no fixed ordering charge asso- 
ciated with any given buy. Various two-period models based 
on different demand probability distributions have been de- 
rived and analyzed, including the dual-period binomial system 
which most closely represents the situation applicable to a 
NARF. In this chapter, several general conclusions are drawn 
from the information previously presented. Additionally, 
recommendations for potential future actions in this area are 
made . 

B. CONCLUSIONS 

The principal results obtained from considering the ap- 
plication of a dual-period binomial inventory model to a NARF 
are listed as follows: 

1. The model identifies a critical inventory quantity 
which can be compared to the previous period closing 
balance to determine both whether additional stock is 
needed and the amount of the stock order if replenish- 
ment action is necessary. 
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2 . The optimal inventory quantity obtained from the two- 
period model is always greater than or equal to the y* 
value determined from a one-period system. As a re- 
sult, there are expected cost benefits to be derived 
from using a two-period approach. 

3. The two-period situation is more difficult to evaluate 
than the single period since expected optimal results 
from the second period must be used in conjunction with 
anticipated first-period demand to determine y*. 

4. The difference between the expected optimal TVC results 
obtained from one and two-period inventory models is 
relatively small. This difference is further reduced 
by the consideration that the production schedule fore- 
cast will generally be less accurate for the second 
period than it will be for the first. 

5. A fairly accurate estimate of the stockout cost rate is 
an essential element in the dual-period formulation. 
This is especially true for high cost items where the 
unit price will be relatively close to realistic stock- 
out cost values. On the other hand, the results ob- 
tained from the dual-period binomial model are not 
particularly sensitive to changes in the surplus cost 
rate . 



C . RECOMMEND AT IONS 

Several recommendations for future actions can be identi- 
fied as a result of the conclusions made in this thesis. 
First, the sensitivity of model results to changes in the 
stockout cost rate indicates that the NARF should attempt to 
quantify all costs incurred as a result of work stoppages so 
that an accurate estimate of P can be made. This is an es- 
sential requirement for identifying the optimal quantity of 
a repair part to stock at a particular point in time, and is 
necessary whether a single or dual-period inventory model is 
used in the RSS. 
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The second recommendation is derived from the conclusion 
that the single-period model provides results which are 
nearly equivalent to those obtained from the dual-period 
system in terras of the expected TVC incurred while using an 
objective function which is much easier to evaluate. Addi- 
tionally, the fact that production schedule forecasts for the 
second period are probably less accurate than for the first 
increases the appeal of the single-period approach. Accord- 
ingly, it is recommended that the single-period binomial in- 
ventory model as discussed in Chapter II of this thesis be 
used by the supporting NSC to determine RSS stock quantities 
for the NARF instead of the dual-period model, pending re- 
ceipt of the results of additional analysis as recommended in 
the following paragraph. 

Finally, the availability of workload forecasts beyond 
the second period must be considered. Although the cost ad- 
vantages obtained from the two-period model over the single- 
period system are small, the possibility exists that these 
benefits are increased by incorporating this additional pro- 
duction information into the model. Therefore, it is recom- 
mended that the two-period model be expanded to consider 
projected rework data spanning up to a procurement leadtime 
in length, and that the optimal results from such a model be 
compared to the single-period results. 
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